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ABSTRACT 



Context. An important open question is the relation between intracluster light and the halos of central galaxies in galaxy clusters. 
Aims. Here we report results from an on going project with the aim to characterize the dynamical state in the core of the Hydra I 
(Abell 1060) cluster around NGC 3311. 

Methods. We analyze deep long-slit absorption line spectra reaching out to ~ 25 kpc in the halo of NGC 331 1. 
Results. We find a very steep increase in the velocity dispersion profile from a central ao = 150 km s -1 to a out — 450 km s -1 at 
R ~ 12 kpc. Farther out, to ~ 25 kpc, a appears to be constant at this value, which is ~ 60% of the velocity dispersion of the Hydra 
I galaxies. With its dynamically hot halo kinematics, NGC 331 1 is unlike other normal early-type galaxies. 

Conclusions. These results and the large amount of dark matter inferred from X-rays around NGC 331 1 suggest that the stellar halo 
of this galaxy is dominated by the central intracluster stars of the cluster, and that the transition from predominantly galaxy-bound 
stars to cluster stars occurs in the radial range 4 to 12 kpc from the center of NGC 3311. We comment on the wide range of halo 
kinematics observed in cluster central galaxies, depending on the evolutionary state of their host clusters. 
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1. Introduction 

An important open question is the physical and evolutionary 
relation between the intracluster light (ICL) and the extended 
halo of the brightest cluster galaxies (BCGs), whether they are 
truly independent components or whether the former is a radial 
extension of the la tter. Using a sample of 683 SDSS clusters 
IZibetti et alJ (I2005I) found a surface brightness excess with re- 
spect to an inner R 1 / 4 profile that characterizes the mean profile 
of the BCGs, but it is not yet known whether this cD envelope is 
simply the central part of the cluster's diffuse light component, 
or whether it is distinc t from the ICL and part of the host galaxy 
dGonzalez et alJl2005l) . 

In the Southern hemisphere, the cD galaxy NGC 3311 and 
the giant elliptical NGC 3309 dominate the central region of the 
Hydra I cluster, an X-ray bright, non-cooling flow, medium com- 
pact cluster with a v elocity dispersion (THydraX — 784 kms -1 
(Mis geld et alJl2008h . The X-ray observations show that the hot 
intracluster medium centered on NGC 331 1 has a fairly uniform 
distribution of temperature and metal abundance from a few kpc 
out to a radius of 230 kpc (iTamura et~al1l2000l : Yamasaki et alJ 
l2002tlHavakawa et al.ll2004lT 2006). Given the overall regular X- 
ray emission and temperature profile, the Hydra I cluster is con- 
sidered as a prototype of an evolved and dynamically relaxed 
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cluster; it is therefore a suitable candidate for a dynamical study 
of a relaxed extended stellar halo around a BCG 

The primary goal of this work is to establish the dynamical 
state of the stellar halo of NGC 3311. We use long-slit spectra 
to uncover the kinematics in the halo region of NGC 3311 out to 
~ 25 kpc from its center. In Sect. [2] we present observations with 
FORS2 at VLT and the GEMINI GMOS archive data, which 
we reanalyze. We describe the data reduction and the kinematic 
measurements in Sect. [3] The newly measured halo kinematics 
and their implications are discussed in Sect. |U and our conclu- 
sions are summarized in Sect. [5] 

We adopt a distance to NGC 331 1 of 51 Mpc (NED), equiv- 
alent to a distance modulus of 33.54 mag. Then 1" corresponds 
to 0.247 kpc. 



2. Observations and archive data 

VLT FORS2 long slit spectra - The long-slit spectra were ob- 
tained during the nights of 2009 March 25-28, with FORS2 on 
VLT-UT1. The instrumental setup had a long-slit l'/6 wide and 
6.'8 long, Grism 1400V+18, with instrumental dispersion of 0.64 
A pixel -1 and spectral resolution a = 90 km s -1 , and a spatial 
resolution along the slit of 0'.'252 pixel -1 . The seeing during ob- 
servations ranged from 0'.'7 to l'/2. The wavelength coverage of 
the spectra is from 4655 A to 5965 A, including absorption lines 
from H^, Mgl (AA5167, 5173, 5184 A) and Fe I (AA5270, 5328 
A). We obtained eight spectra of 1800 sec each, for a total ex- 
posure time of 4 hrs. In the FORS2 observations, the long slit is 
centered on the dwarf galaxy HCC 26 at a = 10h36m45.85s 
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Fig. 1. Optical DSS T x 7' image centered on NGC 331 1 in the 
Hydra I cluster. The relative positions of the GMOS slit (5.'5, 
blue line) and FORS2 slit (6.'8, red line) are illustrated. Green 
and yellow sections on the FORS2 slit indicate regions where 
average spectra are extracted. The adjacent numbers specify the 
radial distances of their light-weighted mean positions from the 
center of NGC 3311. The center of the FORS2 slit coincides 
with the position of the dwarf galaxy HCC 26 and is marked by 
a black circle. North is up and East to the left. 

and 5 = -27d31m24.2s (J2000), with a position angle of 
P.A.=142°; HCC 26 is seen in projection onto the NGC 3311 
halo. The position of the FORS2 long slit is shown in Fig.Q] Its 
center is located at P.A.=64° with respect to NGC 331 1, approx- 
imately along the major axis of the galaxy. 

Gemini GMOS-South long slit spectra - We use Gemini 
archive long-slit spectra in the wavelength range from 3675 A to 
6266 A observed with the B600 grating, a dispersion of 0.914 
A pixel -1 , a spectral resolution of a = 135 kms -1 , and a 
spatial scale of 0'.'146 pixel -1 ; a detailed description of the in- 
strumental setup is presented in Loubser et al. (2008). The see- 
ing was typically in the range from 0'.'6 to l'/2. We target the 
same absorption lines as for the FORS2 spectra, i.e. H^j, Mgl 
(AA5167, 5173, 5184 A) and Fe I (AA5270, 5328 A). The 0'.'5 
wide and 5.'5 long Gemini slit is centered on NGC 3311, at 
a = 10h36m42.74s and S = -27d31m41.3s (J2000), along 
RA. = 63°, the direction of the galaxy major axis. Its position is 
shown in Fig.Q] 

3. Data reductions 

The data reduction of the FORS2 long slit spectra is carried out 
in I RAF. After the standard operations of bias subtraction and 
flat-fielding, the spectra are registered, co-added and wavelength 
calibrated. The edges of the FORS2 slit reach well into sky re- 
gions, which are then used to interpolate the sky emission in the 
regions covered by the stellar spectra. 

In the low surface brightness regions, spectra are summed 
along the spatial direction in order to produce one-dimensional 
spectra with an adequate S/N ratio (> 20 per A). Four indepen- 
dent one-dimensional spectra are extracted along the slit where 
the light is dominated by the halo of NGC 3311; of these, two 
are from regions north and two from regions south of HCC 26, 
respectively. We extract spectra from slit regions of ~ 31" x l'/6 
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Fig. 2. Kinematic fits with PPXF of the spectra extracted at —47" 
(VLT FORS2) and at -34" (Gemini GMOS). In black we dis- 
play the galaxy spectra, in green the wavelength range excluded 
from the fit because of sky residuals, and in red the best-fit- 
broadened template model. All spectra are normalized to the 
value of the best-fit model at 5100 A. 



and ~ 25" x 1"6 at distances of about 54" and 47" from the 
center of NGC 3311, and of - 58" x 1"6 and - 36" x 1"6 at 
central distances of about 100" and 67". We properly mask the 
spectra of foreground stars in those areas. 

The data reduction for the GMOS long slit spectra is car- 
ried out independently here, also in I RAF and with the stan- 
d ard tasks in the Gem ini package. The procedure is described 
m lLoubseret all ([20081 for the wavelength calibration and back- 
ground subtraction; also in this case the edges of the slit are used 
to interpolate the sky emission in the regions covered by the stel- 
lar continuum. Because our goal is to sample the kinematics well 
into the halo, the one-dimensional spectra for the absorption line 
measurements are summed along the slit direction so that a min- 
imum S/N ~ 20 per A is obtained in each radial bin, out to a 
radial distance of about 40" from the center of NGC 3311. 

Stellar kinematics - The stellar kinematics is measured from 
the extracted ID spectra in the wavelength range 4800 A< A < 
5800 A, using both the "penal ized pixel-fitting" method (PPXF, 
ICappellari & Emsel lem (2004)) and the Fourier correlation quo- 
tient (FCQ) method dBenderll 19901) . in order to account for pos- 
sible systematic errors and template mismatch. 

In th e PPXF method, stellar temp late stars from the MILES 
library dSanchez-Blazquez et alJ 12007) are combined to fit the 
one-dimensional extracted spectra; the rotational vel ocity, the 
veloci ty dispersion and Gauss-Hermite moments (e.g. [Gerhard] 
(1 19931) ) are measured simultaneously. Fig. [2] shows two of the 
extracted spectra and the broadened templates fit by PPXF. In the 
FCQ method, the rotational velocity and velocity dispersion are 
derived for each extracted one-dimensional spectrum by assum- 
ing that the LOSVD is described by a Gaussian plus third- and 
fourth- order Gauss-Hermite functions. Before to the fitting pro- 
cedure the MILES template spectra are smoothed to the GMOS 
and FORS2 spectral resolution with the measured broadening 
offsets. While FCQ provides error estimates along with the mea- 
sured kinematics, errors for the PPXF kinematic parameters are 
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Fig. 3. Major axis line-of-sight velocity and velocity dispersion 
profiles for NGC 3311 (P. A. = 6 3°). The open light b lue tri- 
angles are the values published by lLoubser et al.l d2008l) . based 
on Gemini-South GMOS data. The black asterisks are our inde- 
pendent measurements from these GMOS (archival) spectra, and 
the red asterisks show measurements from the new VLT-FORS2 
spectra. These are weighted averages of three independent mea- 
surements which are obtained with PPXF and FCQ as described 
in Sec. [3] and shown separately as the gray, magenta and green 
diamonds. The FORS2 data points are obtained from averages 
over ~ 25" and ~ 31" in the inner regions and over ~ 36" and 
~ 58" in the outer regions of the slit; see Fig.Q] These off-axis 
measurements are plotted at their light-weighted average radii, 
corrected for projection onto the major axis of NGC 331 1 with 
an isophotal flattening of 0.89. Positive distances are south-west 
from the center of NGC 3311 and negative values are north-east, 
along PA.=63°. 



calculated with a series of Monte Carlo simulations adopting the 
appropriate S/N for each bin. 

Because the stellar populations in cD halos may have differ- 
ent metal abundances and ages from those of the inner regions 

■ 11 1| — 1 ° 

(ICoccato et al.H 2010a b). systematic effects caused by template 
mismatch must be evaluated and accounted for. We therefore ex- 
tract kinematic measurements with PPFX and FCQ as follows: 

1. fit with PPFX the best stellar template from the MILES li- 
brary in the central regions with the highest S/N, and extract 
v and a at all radii, using the same stellar template; 

2. simultaneously fit the best stellar template, v and a in each 
radial bin with PPFX; 

3. adopt the respective best-fit PPXF stellar template to derive 
the LOS VD with FCQ for all radial bins; 

4. finally, average rotational velocities v and velocity disper- 
sions a are computed as weighted means of the three val- 
ues extracted in each radial bin as detailed above. Errors 
for these weighted average values are computed from those 
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measurements, but if the reduced \ = 
is greater than one, they are increased 



by y x 2 in order to take into account systematic differences. 
I.e., e| = 1 1 / t i x X 2 where e^, e s are the errors on the 
individual measurements xi and the weighted mean x, re- 
spectively. 

Mean velocities and velocity dispersions in all radial bins 
are listed in Table [T] and the profiles are sho wn in Fig. |3] 
togeth er with the previous measurements from lLoubser et al.l 
(2008). Table 1, which is available in electronic form, contains 
the following information: source of data (Col. 1), distance from 
galaxy center (Col. 2), PA. (Col. 3), v,a with errors for each 
of the procedures 1.-4. described in the text, in Col. (4-5), (6-7), 
(8-9), and (10-11), respectively. Heliocentric and relativistic cor- 
rections have been applied to the mean velocities. The systemic 
velocity is 3800 km s _1 and has been subtracted. 

In the central region of NGC 3311, our ne w velocity dis- 
persio n profile marginally agrees with that of lLoubser et al.l 
(2008). The new mean line-of-sight velocity measurements 
agree with the syste mic velocity of NGC 3311 obtained by 
iMisgeld et al.1 (120081). but ha ve a systematic offset from the v 
data of lLoubser et al.l ([2008 ), by about 91 km s _1 . The agree- 
ment between the new FORS2 measurements at —47" and the 
revised value at —34" from archive GMOS data gives us con- 
fidence that the systematic effects from wavelength calibration 
offsets, template mismatch, etc., are sufficiently small in the new, 
independent data reductions. However, several tests have con- 
vinced us that the data do not allow us to reliably determine full 
line-of-sight distributions (e.g., /13, hn), which could be used to 
test for subcomponents, which one would expect in particular at 
radii ~ 30" -40". 



4. The kinematics of the NGC 331 1 stellar halo 



The combined new velocity dispersion profile for NGC 3311 
reaches to R rnj = 39" ~ 10 kpc from the center of NGC 331 1 
along the galaxy's major axis (PA. =63°), and to an off-axis dis- 
tance of R = 100" ~ 25 kpc along the FORS2 slit. It shows 
a very unusual steep rise with increasing radial distance from 
the galaxy center: from a central value oo = 150 km s -1 , to 
a = 231 km s _1 at R — 15" ~ 3.7 kpc, and then on to 
a flat u out ~ 450 km s -1 outside R = 47" = 12 kpc. The 
steep outward gradient is supported by two i ndependent data sets 
and data reductions. The measurements of Loubser et al. (2008) 
near the galaxy center had already hinted at a positive gradient 
from 190 km s _1 at R = 5" to ~ 24 km s" 1 at a rad ius of 
R = 10", and data shown in Fig. 1 of lHauetall d2004 reach 
~ 300 km s _1 at ~ 25". With the new data we now have very 
clear evidence of a dynamically hot stellar halo in NGC 331 1. 

To put the extremely rapid rise of the velocity dispersion pro- 
file of NGC 33 1 1 in context, we compare its kinematic properties 
with those of earl y-type galaxy (ETG) halos mapped using plan- 
etary nebula data dCoccato et al.ll2009l) and with the halos of the 
two Coma BCG galaxies N GC 4889, NGC 4874 from deep ab- 
sorption line spectroscopy (ICoccato et al.l l2010a). Fig. |4] shows 
the mean < V/a >, X-ray luminosity, and total absolute B-band 
magnitude for these galaxies versus their outermost halo velocity 
dispersion. For NGC 331 1, we use a bolometric X-ray luminos- 
ity within 50" ~ 12 kpc, L x = 2.75 x 10 40 erg s" 1 (based 
on the flux in the 0.5-4.5 keV energy range from Yamasaki et al.l 
(200 2J) and corrected to bolometric Lx according to Table 1 of 
lO' Sullivan et al.l (l200lh \ and the total extinction correct ed B- 
band magnitude (12.22) from Ide Vaucouleurs et al.l d 199 lb . For 
the velocity dispersion of NGC 3311, we use the values at the 
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Fig. 4. Properties of the stellar halo of NGC 3311 compared with 
other early-type galaxy halos: mean < V/a > (upper panel), 
total X-ray luminosity (central panel), and B-band total magni- 
tude (lower panel) against stellar velocity dispersion a. Solid di- 
amond, circle, and square show the measured a of NGC 331 1 at 
the center, 15" (~ 3.7 kpc), and 47" (~ 12 kpc). Crosses show 
outermost velocity dispersions from ICoccato et al.l d2009|). an d 
open diamonds for NGC 4889/4874 from lCoccato et al.l(l2010ab . 



center, at 15" (~ 3.7 kpc) and at 47" (~ 12 kpc). Only the cen- 
tral do puts NGC 3311 in the middle of the ETG distribution; 
cr(47") deviates strongly, with a much larger a than expected 
for the given Lx, Bt- 

The natural interpretation for these results is that the outer 
stellar halo of NGC 33 1 1 is dominated by the central intracluster 
star component of the Hydra cluster. This is supported by sev- 
eral pieces of evidence: (i) The steep rise of the er-profile; more 
isolated ETGs all ha ve slightly or even steep falling er-profiles 
dCoccato et al.ll2009l) . (ii) The saturation of a at ~ 12 kpc, out- 
side of which the dynamically hot component dominates com- 
pletely. ct(47") is ~ 60% of the galaxy velocity dispersion in 
the cluster core, (iii) The large amount of dark matter inferred 
from X -ray observations aroun d NGC 331 1 (~ 10 12 M Q within 
20 kpc. lHavakawa et al.1 d2004l) ). 

In rece nt cosmological hyd rodynamic simulations of cluster 
formation, iDolag et ail (1201 Oh applied a kinematic decomposi- 
tion to the stellar particles around cD galaxies. With a double 
Maxwellian fit to the velocity histogram of star particles cen- 
tered on a simulated cD, they were able to separate an inner, 
colder Maxwellian distribution associated with the central gal- 
axy, and an outer, hotter component of stars that orbit in the clus- 
ter potential. For both components they derived radial density 
profiles and, fitting Sersic profiles, found that the inner stellar 
component is much steeper than the outer diffuse stellar com- 
ponent. A comparison with these simulations indicates that the 



steep velocity dispersion gradient in the halo of NGC 3311 traces 
the transition from central galaxy stars to the diffuse intracluster 
stellar component. In the NGC 331 1 halo, the transition between 
the two occurs at smaller radii than in other BCGs in nearby clus- 
ters, in the range between 4 and 12 kpc. 

NGC 3311 appears to have a si milar halo as the c D galaxy 
NGC 6166 in the Abell 2199 cluster dKelson et al.ll2002l) . whose 
CT-profile rises to cluster values at R ~ 60 kpc. But NGC 3311 
is even more extreme; it is a fairly small galaxy, based on its 
central cto = 150 kms -1 , and it is already dominated by the 
surrounding cluster component at R ~ 12 kpc. Presumably, this 
is because the core of the "relaxed" Hydra cluster has had time 
to collapse onto the galaxy. For comparison, the two BCG gal- 
axies in t he Coma cluster core , which have a nearly constant 
CT-profile dCoccato et al.ll2010 a). may be in the middle of an on- 
going merger (Ger hard et al.ll2007h . so that their previous sub- 
cluster halos would have been stripped and a new cluster halo 
could have been built only after the merger was completed; and 
in the outer halo of t he more isolated M8 7, the velocity disper- 
sion appears to drop (iDohertv et al.ll2009b towards the edge. 



5. Conclusions 

Based on two independent long-slit data sets and reductions, we 
find a steep gradient in the velocity dispersion profile of the cen- 
tral galaxy NGC 3311 in the Hydra I cluster, from cto — 150 
km s _1 to a ou t — 450 km s _1 outside 12 kpc (60% of the ve- 
locity dispersion of the galaxies in the surrounding cluster). 

The new data provide evidence that NGC 3311 is a fairly 
small galaxy dominated by a large envelope of intracluster stars 
already beyond R ~ 12 kpc, whose orbits are dominated by the 
cluster dark matter potential. Comparison with other BCG gal- 
axies shows a wide range of dynamical behavior in their halos. 
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Table 1. Measured mean velocities and velocity dispersions for NGC 3311. For details see text. The galaxy's systemic velocity V sys = 3800 km s 1 has been obtained by a 
linear fit to the velocities in the central 20" and has then been subtracted from the measurements. This value includes heliocentric and relativistic corrections. 
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67 ± 13 


456 ± 


15 


87 ± 13 


447 ± 


15 


124 ±8 


467 ± 14 


105 ± 17 


457 ±8 


FORS2 




-67 


6° 


87 ±53 


471 ± 


49 


66 ±53 


454 ± 


49 


97 ±55 


476 ± 97 


83 ±31 


464 ± 33 


FORS2 




-100 


114° 


14 ±21 


437 ± 


35 


47 ±21 


432 ± 


35 


50 ±29 


405 ± 51 


35 ± 13 


429 ± 22 



Notes - Col. 1: instrument. Col. 2: Radial distance from center of NGC 3311. Col. 3: Position angle of data with respect to NGC 3311's center. Col. 4: Velocity measured with PPXF (using 
template determined at R = 0"), relative to the galaxy systemic velocity. Col. 5: Velocity dispersion measured with PPXF (using template determined at R = 0"). Col. 6: Velocity measured 
with PPXF (free template), relative to the galaxy systemic velocity. Col. 7: Velocity dispersion measured with PPXF (free template). Col. 8: Velocity measured with FCQ, relative to the galaxy 
systemic velocity. Col. 9: Velocity dispersion measured with FCQ. Col. 10: Weighted average velocity. Col. 11: Weighted average velocity dispersion. 



